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Abstract. The interacting boson model (IBM) of Arima and lachello is applied to calculate nuclear 
charge radii and electric monopole transitions of even-even nuclei in the rare-earth region. Consis- 
tent operators are used for the two observables. A relation between summed Ml strength and p (EO) 
values is pointed out. 
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1. INTRODUCTION 

According to the geometric model of Bohr and Mottelson [1], a nucleus with an elip- 
soidal equilibrium shape may undergo oscillations of two different types, /3 and y. While 
7 vibrations are an acknowledged feature of deformed nuclei, such is not the case for the 
/3 kind. A careful analysis of the observed properties of excited 0+ states seems to indi- 
cate that very few indeed satisfy all criteria proper to a /3 -vibrational state [2] . In particu- 
lar, very few decay to the ground state by way of an EO transition of sizable strength [3], 
as should be for a /3 vibration [4]. It is therefore not surprising that alternative inter- 
pretations of excited 0+ states in deformed nuclei, either as pairing isomers (see, e.g., 
Refs. [5, 6]) or through shape coexistence and configuration mixing (see, e.g., Ref. [7]) 
have gained advocates over recent years. 

In this contribution it is shown that a collective interpretation of the ground and 
excited 0+ states with the interacting boson model (IBM) of Arima and lachello [8] 
can account for the data on charge radii and EO transitions, as observed in the rare-earth 
region. Full details of this work can be found in Refs. [9, 10]. In addition, a correlation 
with summed Ml strength is pointed out. 



2. CHARGE RADII AND ELECTRIC MONOPOLE TRANSITIONS 

The probability for an EO transition to occur between an initial state \i) and a final state 
I/) can be written as the product of an electronic factor Q. and the square of a nuclear 
factor p (EO) , the latter being equal to 



with R = roA^/^. In first approximation the EO matrix element equals 



{f\fp{EO)\i) = {f\e'£rl\i). (2) 

where the sum is over all protons in the nucleus [11]. In the nuclei considered here not 
all protons can be treated explicitly and recourse should be taken to effective charges en 
and for the neutrons (n) and protons (p), leading to the generalized expression [12] 




On the other hand, the mean-square charge radius of a state \s) is given by 

(.|7p(r2)|.) = l(.|j:r2|5). (4) 

This is an appropriate expression insofar that a realistic A-body wave function is used 
for the state I*). This is often impossible and an effective charge radius operator r(r^) 
should then be taken. The generalization of the expression (4), similar to the one carried 
out for the EO operator, can therefore be written as 

(r2), ^ {s\f{r')\s) = \ {s\ £ e,rl\s). (5) 

CniV -r cpZ, 

The basic hypothesis of the present study is to assume that the effective nucleon charges 
in the charge radius and EO transition operators are the same. If this is so, comparison 
of Eqs. (3) and (5) leads to the operator relation 

f{m) = {enN^e^Z)f{r^). (6) 

This is a general relation between the effective operators to be used for the calculation of 
charge radii and EO transitions. Equation (6) can, in principle, be tested in the framework 
of any model and here the implied correlation is tested in the framework of the IBM [8]. 
This requires that all states involved [i.e., \i) and |/) in Eq. (1), and I*) in Eq. (5)] are 
collective in character and can be described by the IBM. 

In the IBM-1 the charge radius operator is taken as the most general scalar expression, 
linear in the generators of U(6) [13], 

f(r2) = (r2), + aArb + 77^, (7) 

where A'b is the total boson number, fid is the J-boson number operator, and a and r\ 
are parameters with units of length^. The first term in Eq. (7), (r^)c, is the square of 
the charge radius of the core nucleus. The second term accounts for the (locally linear) 
increase in the charge radius due to the addition of two nucleons {i.e., neutrons since 



isotope shifts are considered in this study). The third term in Eq. (7) stands for the 
contribution to the charge radius due to deformation. It is identical to the one given 
in Ref. [13] but for the factor 1/A^b- This factor is included here because it is the 
fraction {nj) /Nh which is a measure of the quadrupole deformation (jSl in the geometric 
collective model) rather than the matrix element (hj) itself. 

Two quantities can be derived from charge radii, namely isotope and isomer shifts. 
The former measure the difference in charge radius of neighboring isotopes. For the 
difference between even-even isotopes one finds from Eq. (7) 

A{.^) . (^><f > - (^><|) = i«i + n - <|>ir') ■ («) 

Isomer shifts are a measure of the difference in charge radius between an excited (here 
the 2^) state and the ground state, and are given by 

5(,^>.(.^)<f-(.^)i|l = i((«4'-(«.>l|'). (9) 

Once the form of the charge radius operator is determined, the EO transition operator 
follows from Eq. (6). In the IBM-1 the EO transition operator is therefore 

f(E0) = (enA^ + epZ)7]^. (10) 

Since for EO transitions the initial and final states are different, neither the constant (r^)c 
nor aNb in Eq. (7) contribute to the transition, and its p (EO) value equals 

p(EO;/^/) = ^^^^^A|(/|^,|,)|, (11) 



3. APPLICATION TO THE RARE-EARTH REGION 



To test the relation between charge radii and EO transitions, proposed in the previous 
section, a systematic study of all even-even isotopic chains from Ce (Z 58) to W 
(Z = 74) is carried out. This analysis requires the knowledge of structural information 
concerning the ground and excited states which is obtained by adjusting an IBM-1 
Hamiltonian to observed spectra in the rare-earth region. The details of the energy 
calculation can be found in Ref. [10]. The procedure is closely related to the one 
followed by Garcia-Ramos et al. [14] and yields a root-mean square deviation for an 
entire isotopic chain which is typically of the order of 100 keV. 

Isotope shifts A(r^), according to Eq. (8), depend on the parameters \a\ and r\ in the 
IBM-1 operator (7). The parameter |a| is adjusted for each isotope series separately, 
while r\ is kept constant for all isotopes, r\ = 0.50 fm^. The resulting isotope shifts 
are shown in Fig. 1. The peaks in the isotope shifts are well reproduced in all isotopic 
chains with the exception of Yb. The largest peaks occur for ^^^^^^''Sm, ^^^^^^^Gd, and 
156-1 54 that is, for the difference in radii between A' = 90 and N — SS isotopes. The 
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FIGURE 1. Isotope shifts A{r'^) in the rare-earth region. The points are the experimental values and the 
Unes are calculated with Eq. (8). References to experiments are given in Ref. [10]. 
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FIGURE 2. Isomer shifts 5 (r2 ) in the rare-earth region. The open symbols are the experimental values 
and the full symbols are calculated with Eq. (9). References to experiments are given in Ref. [10]. 



peak is smaller below Z = 62 for Ce and Nd, and fades away above Z = 66 for Er, Yb, Hf, 
and W. The calculated isotope shifts broadly agree with these observed features but there 
are differences though. Notably, the calculated peak in the Sm isotopes is much broader 
than the observed one, indicating that the spherical-to-deformed transition occurs faster 
in reality than it does in the IBIVI-I calculation. 

A further test of the calculated charge radii is obtained from isomer shifts 5(r^), 
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FIGURE 3. Electric monopole p^(E0;7f Jf) values (times 10^) in the rare-earth region for = 0+ 
and 2+. The open symbols are the experimental values and the full symbols are calculated with Eq. (11). 
References to experiments are given in Ref. [10]. 



depending only on rj [see Eq. (9)]. The isomer shifts known experimentally are shown 
in Fig. 2. The data are more than 30 years old and often discrepant, in particular in 
the Yb and W isotopes. Nevertheless, a clear conclusion can be drawn from the isomer 
shifts measured in the Sm and Gd isotopes: they are an order of magnitude smaller in 
the deformed than they are in the spherical region. 

The question is now whether the value 7] = 0.50 fm^, obtained from isotope and 
isomer shifts, reproduces the EO transitions observed in the rare-earth nuclei. In Fig. 3 
the available EO data for 0^ — > 0^ and 2f 2^ transitions in the rare-earth region 
are compared with the results of the IBM-1 calculation. It is not the intention to give 
a detailed comparison for all transitions (which can be found in Ref. [10]) but the 
overall impression is that the present approach succeeds in reproducing the correct order 
of magnitude for p^(EO). However, some discrepancies can be observed in heavier 
nuclei and especially concern ^^^Yb and 1^2-184-^ ^ possible explanation is that the 
p^(EO) measured for these nuclei is not associated with collective states. This seems 
to be the case in ^^^Yb where several p^(EO) have been measured, none of which is 
large. Only in the W isotopes does it seem that the observed EO strength is an order of 
magnitude smaller than the calculated value. It is known that these nuclei are in a region 
of hexadecapole deformation [15] and this may offer a qualitative explanation of the 
suppression of the EO strength [10]. 

4. RELATION TO SUIMIMED IMl STRENGTH 

After the presentation of this work during the conference, J.N. Ginocchio asked a 
question concerning a possible correlation with summed Ml strength to the scissors 



mode. The argument goes as follows. It is known from the work of Ginocchio [16] that 
the summed Ml strength from the ground state to the scissors mode (for a review on the 
latter, see Ref. [17]) is related to the ground-state matrix element of the <i-boson number 
operator Ud, 

y{Ml-Ot ^ i;) = 4(^v -g.)^^^^^^(0+|n.|0+), (12) 

where n (z) is the number of valence neutron (proton) particles or holes, whichever is 
smaller, and and g^i are the g-factors of the neutron and proton bosons. Note that 
n + z = 2N\,. The relation (12) was used in Ref. [18] to establish a connection between 
the summed M 1 strength and the charge radii of the Nd, Sm, and Dy nuclei. It is therefore 
natural to ask whether in turn a connection exists between the summed Ml strength and 
p (EO) values. 

Such a connection can indeed be established. The EO operator is directly proportional 
to n^, unlike the charge-radius operator which involves additional terms [see Eq. (7)] 
which complicate the relation between the summed Ml strength and charge radii. On 
the other hand, the matrix element of ha appearing in the sum rule (12) is diagonal while 
a p (EO) value involves a non-diagonal matrix element of hd. A relation between the two 
matrix elements can nevertheless be obtained in the SU(3) limit where they are [19] 



,„+, 4iVb(iVb-l) + I 



8(iVb-l)"yVb(2iVb + l) 
9(2A^b-3)(2A^b-l)2 



1/2 



(13) 



From these expressions the ratio of matrix elements can be derived, resulting in the 
following relation, valid in the large-A^b limit: 

5(M1;0+ ^ 1+) ^ ^{gy-gnf'ig{N,Z,n,z)p{EO;0+^0+), (14) 

where g{N, Z, n, z) is the function 

,(N,Z,n,.)=<!^^. (15) 

CnN + epZ ^Jn-\-z 

Note that in the SU(3) limit only one l"*" state is excited and only the j8 -vibrational state 

decays by EO to the ground state, as is indicated in Eq. (14). 

It should be emphasized that the relation (14) is valid only in the SU(3) limit which 
might jeopardize its use in transitional nuclei. One may nevertheless attempt to apply it 
to the entire rare-earth region. The ratio Tq/t; = 3.08 and the effective charges = 0.5e 
and = e have been determined from a fit to radii [10]. The only remaining constant 
in Eq. (14), gy — gn, can be obtained by adjusting the expression (12) to the observed 
summed Ml strength in rare-earth nuclei [20], as shown in Fig. 4. This, incidentally, 
constitutes an additional test of the matrix elements (0^|nrf|0^) as they are obtained 
in the IBM-1 fits to the different isotope series. An overall satisfactory agreement is 
obtained with the value \gy — gn\= 0-83 /^n but it is seen that the calculated Ml strength 
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FIGURE 4. Summed Ml strength Y,fB{Ml;Qi 1^) in the rare-earth region. The open symbols are 
the experimental values [20] and the full symbols are calculated with Eq. (12). 
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FIGURE 5. Correlation between the p(EO;0^ 0^) value and the summed Ml strength 
J^-fB{Ml;0^ 1^) for the eight nuclei in the rare-earth region where both properties are known. The 
function g{N,Z,n,z) is defined in Eq. (15). The slope of the line is 9{gv — gK^rQ/ [Ani]) jl^. 



in the Sm isotopes is well below its observed value, again indicating that the calculated 
transition is too slow for this isotope series. 

The correlation (14) can now be tested and is shown in Fig. 5 for the eight nuclei 
in the rare-earth region where both EO and Ml properties are known (I50,i52,i54gj^^ 
154,156,158q^^ I66gj.^ ^j^^ niyi^y jj^g ^^^Yh point is conspicuously off the line, reinforc- 
ing the earlier remark about EO strength in this nucleus. The ^^^Er point follows from a 
recent experiment [21] where the fourth = 0"*" level at 1934 keV has been identified as 
the band head of the j8 -vibrational band with a sizable EO matrix element to the ground 
state. 



5. SUMMARY 



A consistent description of nuclear charge radii and EO transitions was proposed and the 
vaUdity of this approach was tested with the interacting boson model applied to even- 
even nuclei in the rare-earth region (58 < Z < 74). An additional correlation between 
p{EO) values and the summed Ml strength to the scissors mode, valid in the SU(3) 
limit, was pointed out and tested in the rare-earth region. It would be of interest to 
derive a more general relation between these two properties, applicable to transitional 
and deformed nuclei. 
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